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Abstract 

We show that the [i problem and the strong CP problem can be resolved in the context 
of the gauged U(1)r symmetry, realizing an automatic Peccei-Quinn symmetry. In 
this scheme, right-handed neutrinos can be introduced to explain small Majorana 
or Dirac neutrino mass. The U(1)r D-term mediated SUSY breaking, called the 
U{1)r mediation, gives rise to a specific form of the flavor-conserving superpartner 
masses. For the given solution to the \x problem, electroweak symmetry breaking 
condition requires the superpartners of the Standard Model at low energy to be much 
heavier than the gravitino. Thus dark matter candidate can be either gravitino or 
right-handed sneutrino. In the Majorana neutrino case, only gravitino is a natural 
dark matter candidate. On the other hand, in the Dirac neutrino case, the right- 
handed sneutrino can be also a dark matter candidate as it gets mass only from SUSY 
breaking. We discuss the non-thermal production of our dark matter candidates from 
the late decay of stau and find that the constraints from the Big Bang Nucleosynthesis 
can be evaded for a TeV-scale stau mass. 
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1 Introduction 



N — 1 supersymmetry (SUSY) contains a continuous U(1)r group transforming differ- 
ent supercharges and thus distinguishing between bosonic and fermionic components of 
superfields. A discrete subgroup of it, called R-parity, can remain respected after SUSY 
breaking, and commonly used to explain the stability of proton. In local supersymmery 
(supergravity) , such an R-symmetry, either discrete or continuous, is gauged and should 
respect anomaly-free condition PQ[2]. A gauged R-symmetry controls all the fields in a 
theory and thus its anomaly-free condition is very restrictive. Due to this property, R- 
symmetries can provide a powerful tool for phenomenological applications such as U(1)r 
as a family symmetry [3], as a resolution to the \i problem and B/L conservation jl], as an 
origin of supersymmetry breaking [5] and the D-term inflation [6]. More recently, an U(1)r 
mediated supersymmetry breaking model has been constructed based on a six-dimensional 
flux compactification [7\ and its phenomenological application was investigated by some 
of the authors [H[9]. In 6D compactifications, even if the hidden sector is geometrically 
separated from the visible sector, the hidden sector SUSY breaking is generically not se- 
questered [ID] , as compared to the 5D counterpart. However, it has been shown that the 
moduli F-term contribution to the soft mass is cancelled by the hidden F-term contribution 
such that the U(1)r mediation can be dominant [EJ|9]. 

In this paper, extending the previous studies [9], we propose a resolution of the \i 
problem along the line of Ref . |TT] , realizing also the axion solution to the strong CP prob- 
lem [12]. At the same time, the observed neutrino masses and mixing can be explained 
by introducing three right-handed neutrinos, which can form Majorana neutrinos by the 
usual seesaw mechanism (at the intermediate axion scale of order 10 10 ^ 12 GeV or at the 
TeV scale) or Dirac neutrinos with tiny neutrino Yukawa couplings. In this framework, we 
find that the electro-weak symmetry breaking (EWSB) condition requires a peculiar su- 
perparticle mass spectrum: all the superparticles in the minimal supersymmetric standard 
model (MSSM) sector have masses of the TeV scale whereas the gravitino or right-handed 
sneutrino masses can be around 100 GeV. This is in contrast to the previous study of the 
U(1)r phenomenology [S] where the gaugino masses are assumed to be comparable to the 
scalar soft masses at the GUT scale and the fi and Bfi terms are assumed to be given such 
that the EWSB condition is satisfied. 

As the superpartners in the MSSM turn out to be heavy, either the gravitino or the 
right-handed (RH) sneutrino can be a natural dark matter candidate. We focus on the 
non-thermal productions of dark matter depending on the nature of neutrinos and impose 
the Big Bang Nucleosynthesis (BBN) constraints. When neutrinos are of Majorana type, 
it is typical that gravitino is the lightest supersymmetric particle (LSP) and stau is the 
next LSP (NLSP). In this case, the stau decay after the freezeout becomes a dominant 
source for the non-thermal production of gravitino relic density. On the other hand, when 
neutrino is of Dirac type, RH sneutrino can be the LSP and then gravitino is the NLSP. 
Because of the tiny Dirac neutrino Yukawa couplings, the thermal production of the RH 
sneutrino from the decay of heavy superparticles is usually suppressed. Then, the decay 
of stau can be a dominant source for the RH sneutrino relic density. We point out that 
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in both gravitino and RH sneutrino dark matter scenarios, the BBN problem coming from 
the late decay of stau can be avoided for the TeV-scale stau mass as required from the 
EWSB conditions. 

The paper is organized as follows. First, we present a 4D effective theory with gauged 
.R-symmetry where the hidden sector SUSY breaking is introduced and the visible sector 
contains no additional representations under the standard model gauge group other than 
the MSSM content. In the next section, we show that the soft mass parameters at the 
GUT scale are determined from the U(1)r mediation and the fj, term and the neutrino 
masses are dynamically generated with nonzero singlet vacuum expectation values (VEVs). 
Consequently, deriving the low-energy SUSY spectrum, we discuss the relic density of dark 
matter and the BBN constraints, depending on the nature of neutrinos. Finally, conclusion 
is drawn. There are four appendices comprised of the minimization of the singlet scalar 
potential, the identification of the axion for multiple scalar VEVs, the discussions on 
the effective -R-parity and PQ symmetry violating terms induced after the .R-symmetry 
breakdown. 



2 A SUSY model with gauged R-symmetry 

We first review the consistency of the Fayet-Iliopoulos(FI) term in supergravity and explain 
the 4D effective gauge supergravity recently derived from a six-dimensional flux compact- 
ification. For generation-independent .R-charges and renormalizable Yukawa couplings in 
the MSSM, the gauged U(1)r is anomaly-free up to a Green-Schwarz counter term. We 
obtain the /i term and the Majorana/Dirac neutrino mass terms through the superpoten- 
tial terms with singlets of intermediate-scale VEVs. In both Majorana and Dirac neutrino 
cases, we present the representative models where negative -R-charges are assigned for the 
MSSM scalar partners except for the Higgs doublets. 



2.1 Fayet-Iliopoulos term in 4D supergravity 

It has been recently pointed out that a U(l) gauge theory with constant FI term can 
be consistent with supergravity, provided that an exact global symmetry is present [13]. 
However, any global symmetry is believed to be broken in quantum gravity, leading to 
a conclusion that there is no consistent U(l) with a constant FI term. Here we review 
the FI term in supergravity focusing on the chiral compensator formalism discussed by 
Komargodski and Seiberg in Ref. [13] and comment on the consistency of the FI term for 
the gauged U(1)r. 

The general 4D gauged supergravity action in Weyl compensator formalism is 

S = Jd 4 x [d 4 0E(-3CV^ v «/ 3 C e -*o(* i Vo,# i /o)/3 ) + J d ^SC 3 W(^/C) + h.c] (1) 

where C is the compensator superfield, which becomes C = C + 9 2 F C in super- Weyl gauge, 
E and £ are the full and chiral superpace measures, respectively, Vr is the U{1)r vector 
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superfield, and qr is the U(1)r gauge coupling. For the gravitino of i?-charge +1, the 
constant Fayet-Iliopoulos(FI) term for the gauged U(1)r is quantized as £ = 2. The above 
supergravity action can be made U(1)r gauge invariant by a super- Weyl transformation 
j6j[8]. A construction of the gauged U(1)r invariant action in 4D supergravity has been 
originally done in Ref. [14] . 

The U{1)r transformation, which should not be confused with the gauged R-symmetry, 
is defined such that all the chiral superfields have R-charge 2/3. We note that the fermionic 
superpartners have R-charge —1/3. The original theory with general superpotential W(<&j) 
is made R-symmetric by adding an additional chiral superfield C with R-charge 2/3. 

Now we consider a local U(1)q without FI term under which the £7(1) charges are 
Q[$i] = Qi and Q[C] = 0. After introducing a nonzero FI term for the U(1)q in the 
Kahler potential, we need to make charge shifts, obtaining a new local U(1)q: U(l) charges 
are Q[C] = — £/3 and Q[&i] = Qi + Ri£/2 — £/3 where Ri are new R-charges satisfying 

Ri — 2 for chiral superfields appearing in each term of the superpotential. 

After C gets a nonzero vacuum expectation value, C = = Mp, the U(1)q and the 
U(1)r is broken down to a gauged R-symmetry, U(1)q + |?7(l)jj = U(1) R . The R-charges 
of this gauged R-symmetry are R[C] = 0, R[Q>i] = Qi + Ri^/2 = Ri. The R-charges of the 
fermionic superpartners of $j are R[ipi] = Qi + (Ri — 1)6/2. Therefore, since £\ Ri — £ and 
Y2i Ri = 2 for chiral superfields appearing in each term of the superpotential, one can draw 
a conclusion that ^ Qi = J2i ^» ~~ (6/2) J2i Ri = 0> so there appears a global symmetry 
for the general superpotential. 

This result is due to the assumption that there is a local U(1)q in the limit of a vanishing 
FI term. In order to construct the U(l) theory with nonzero FI term, however, one only 
has to start with U(1)q symmetry instead of local U(1)q. Furthermore, when the FI term 
is quantized as required for charge quantization, there is no limit of a vanishing FI term. It 
has been shown that a consistent 4D supersymmetric vacuum with gauged R-symmetry can 
be obtained below the compactification scale in six- dimensional gauged supergravity [TIE] - 
In this case, the quantization of the FI term is originated from the flux quantization in 
extra dimensions. For instance, when £ = 2, we only have to start with the U(1)q having 
charges Qi = Ri — | in terms of the R-charges of the gauged R-symmetry, while the U(1)q 
symmetry does not need to be imposed. The charges of the would-be global symmetry 
U(1)q are Qi = Ri — Ri, so U(1)q would be unbroken only if one can find Ri's satisfying 

Ri = 2 at all orders. However, it is also possible to have £\ R j = 2 but Ri ^ 2 for 
higher dimensional terms in the superpotential. 

As will be discussed in later sections, the PQ symmetry appears in our model and 
it is nothing but a global symmetry with new R-charges given by fj = r« + with g; 
being the PQ charges and being the R-charges of the gauged U(1) R . Here, = 2 is 

guaranteed at low orders in the superpotential by the fact that £\ = and J2i r i = 2- 
However, at higher orders, even for = 2, we found that £V $ = or J2i ?i = 2 is 

not satisfied any more. For instance, for X, Y singlets considered in our previous paper, 
the PQ symmetry is an accidental symmetry which holds for the quark/lepton Yukawa 
couplings and the /i term. However, the PQ symmetry is broken by the other Planck-scale 
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suppressed [7(1 ^-invariant interactions. 

2.2 4D effective supergravity from a 6D flux compact ificat ion 

In this section, we consider a concrete form of the Kahler potential and the superpotential 
derived from a flux compactification in six dimensions |8, 9j. The bulk theory is based on a 
6D chiral gauged supergravity constructed by Nishino and Sezgin [15] . The gauged U(1)r 
appears as a partial gauging of the bulk i?-symmetry. In this flux compactification, there 
is a non- vanishing gauge flux along the U(1)r, making the 4D Minkowski space flat while 
the 2D extra dimensions are compactified on the sphere with a wedge cut out 0[16] . There 
are two 3-branes with nonzero equal tension at the poles of the wedged sphere so visible 
sector fields are located at one pole and the hidden sector fields are located at the other 
pole. In order to stabilize the remaining modulus, some bulk dynamics should be taken 
into account too. 

In the 4D effective supergravity, the FI term is given by £ = 2 and part of the Kahler 
potential without FI term is 

K = -ln(~(S + £t)) 

- In (i(T + - 8g R V R ) - Qle^^Qi - Q'^^Q' - ^e^^) 
+M^e~ 2rM9RVR M. (2) 

Here S, T are the moduli that mix the dilaton and the volume modulus, Qi are visible 
brane fields, Q', tp(M) are hidden sector fields living on the hidden brane (in bulk). In our 
model, the gauged R-symmetry is spontaneously broken by a flux compactification in six 
dimensions [8]j9]. So it is non-linearly realized by the axion of the bulk T- modulus and the 
mass of the U(1)r gauge boson is of order the compactification scale. This is manifest in 
the above form of the Kahler potential. 

The superpotential for the modulus and the hidden sector is 

Moduli = Wo + fQ' + -^-e- bS + AV P M 2 + (3) 

where the i?-charges are rqi = 2, r M = —~ and r v = | = ~^r^- Here we introduced 
the uplifting sector parametrized by / and the bulk sector responsible for the gaugino 
condensate, in an explicitly U(1)r invariant fashion. When a Green-Schwarz coupling to 
the T-modulus is responsible for cancelling the i?-symmetry anomalies, the superpotential 
for the gaugino condensate would get a T-dependent factor with S being replaced by S+eT. 
But we assume e«l such that the T-modulus dependence gives a negligible effect on the 
soft scalar masses. 

On the other hand, Wo stands for a nonzero VEV of the superpotential obtained after 
a spontaneous breaking of the U(1)r in the hidden sector. In generalized O'Raifeartaigh 
model with renormalizable interactions, independent of the i?-symmetry breaking, the 
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superpotential VEV is undetermined at tree level as it is proportional to the pseudo- 
moduli [17]. However, it is possible to stabilize the pseudo- moduli at a nonzero value from 
the Coleman- Weinberg potential at one-loop [18]. In this case, the superpotential can get 
a nonzero VEV, Wo ^ 0. However, the .R-symmetry breaking sector giving rise to Wq ^ 
would necessarily break SUSY and generate a positive vacuum energy, because of the 
consistency condition, 2Wq = Y2i r i ( l ) i7^~- For instance, a bulk .R-symmetry breaking field 
$ with R-charge +2 leads to the scalar potential, V$ = ( Re g)( R eT-QtQ) I 1 2 • Because of the 
bound [H] on |W |: 2|W | < f r F with f r = ^rf |0 4 | 2 and F 2 = £\ lfj| 2 , the positive 
vacuum energy could not be cancelled by the perturbative contribution proportional to 
| Vl^o 1 2 for f r ^ Mp. However, since the R-symmetry anomalies are cancelled by a Green- 
Schwarz mechanism, a non-perturbative correction may break the R-symmetry dynamically 
so that we may avoid the no-go theorem based on the perturbative generation of the 
superpotential. In this case, the order parameter of the R-symmetry breaking is now the 
superpotential vev itself, not the vev of a fundamental scalar field. 

The moduli stabilization with the effective superpotential ((3]) has been discussed in 
Ref. [9]. The real part of the T modulus was shown to be stabilized at t ~ 1 mainly by the 
U(1)r D-term. This is due to the cancellation between the constant Fayet-Iliopoulos(FI) 
term present in 4D gauged supergravity and the field-dependent FI term coming from the 
internal gauge flux. On the other hand, the S modulus and M, (f are stabilized by the 
F-terms at the perturbative regime. Q' is also stabilized radiatively due to the supersym- 
metric couplings to heavy fields without introducing additional SUSY breaking sources. 
Then, the resulting F-terms for S, M, ip are negligible while the F-terms for T, Q' are 

F T ~ 2m3/2, F®' ~ y/2m 3 / 2 , respectively, and the t/(l)_R D-term is Dr ~ — . Conse- 

9R i—i 

quently, because of the cancellation between the moduli and hidden-brane F-termaj, scalar 
soft masses are determined dominantly by the U(1)r D-term as 

m 2 ~ rig R D R ~ -rp^ . (4) 

The nonzero soft masses for brane scalars proportional to the -R-charges can be also derived 
directly from the 6D action for the non-supersymmetric flat brane solution with a small 
warping [22]. The warping induced by unequal 3-brane tensions makes the brane- localized 
and flux- induced masses uncancelled [S]. Therefore, as a small U{1)r D-term is generated 
below the compactification scale (or the 4D GUT scale), it gives rise to a visible effect on 
the low-energy phenomenology by contributing to the initial SUSY spectrum at the GUT 
scale, unlike the conclusion of Castano et al in Ref. [2] where it was assumed that the 
.R-symmetry is broken at the Planck scale and the U(1)r D-term vanishes. 

2.3 The U{1)r anomalies 

Assuming that the renormalizable Yukawa couplings are allowed for quarks and leptons in 
the MSSM, the anomaly cancellation conditions for the U{1)r determine the -R-charges of 

1 Due to the absence of sequestering in 6D compactifications [10] . the contact term between hidden 
sector field Q 1 and the visible sector superfields makes the cancellation happen. 
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sfermions in terms of the squark doublet i?-charge q [9] as follows, 

r 16 _ 3_ 26 _ 17 _ 18 

I = —3q , e— — q , u— — q H , 

« 3 > 7 y 21' 7 y 7' 

~ 31 46 ~ 24 60 ~ 24 4 

« = -y?-y, /id = yg+y, K = -—q- -. (5) 

So, there is one parameter family of solutions to the consistent .R-charges. We note that the 
i?-charge of a fermion differs from the one of scalar superpartner by one unit, as I = I — 1. 
We assume that the pure U(1)r anomalies are cancelled by hidden fermion^. On the other 
hand, it has been shown that nonzero {7(1)#-SM mixed anomalies, C a (a = 1,2,3), can be 
cancelled by the variation of a Green-Schwarz term [9]: 

3 1 

£ GS = (ImT)^ a -tr(i^F a ) (6) 

0=1 

where the U(1)r gauge transform of ImT is J^(ImT) = Ag^A^ and k a are related to the 
anomaly coefficients as k a = i 6 ^2 gR with C\ = —15 and C 2 = C 3 = —9. Consequently, 
after a supersymmetric completion of the Green-Schwarz term, the gauge kinetic functions 
for the brane-localized SM gauge fields are modified to 

fa = 4~ + k * T > 00 
9a,0 

where g a are the tree-level SM gauge couplings. For unified tree-level gauge couplings 
with #f o = g\ Q and gf = |^| at the compactification scale, k a = |A; 2 is consistent with 
the favorable choice of sin 2 9w = § at the compactification scale, as there is no exotics 
charged under the SM between the unification scale and the electroweak scale. However, 
given that the R-charges of the matter fields are non-universal in the same GUT multiplet, 
the matter multiplets should appear as split multiplets below the compactification scale as 
in orbifold GUT models. 



2.4 The ii term 

In the presence of the gauged U(1)r, the /i term is forbidden at tree level by the consistent 
i?-charge assignment in (J5J). Therefore, for the resolution of the \x problem, we introduce 
higher dimensional interactions with two singlets, Y,X, in the superpotentiaH: 

W u = —Y 2 H u H d + —YX\ (8) 
" M P M P K J 



2 See the U(1)r anomaly coefficients in the presence of hidden fermions with nonzero R-charges in 
Ref. [9]. There may be an additional R-symmetry breaking in the process of giving hidden fermions 
masses. If this occurs only in the hidden sector, there is no problem with interactions of the additional 
R-breaking fields and the MSSM fields. There might appear also light fermions and one of them could be 
the LSP. 

3 We note that other dimension-5 singlet operators, Y 2 X 2 and Y 3 X, are problematic because there is 
no minimum with nonzero singlet VEV. 
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Here we note that there exists an automatic Peccei-Quinn (PQ) symmetry^ which provides 
the axion solution of the strong CP problem [12]. As shown in Appendix A, in the presence 
of weak-scale soft mass terms, the dimension-5 interaction between singlets gives rise to 
intermediate-scale singlet VEVs so that we can get a weak-scale \i term. Note that this 
realizes the idea of Ref. [TT] by imposing a fundamental U(1)r symmetry. In the same 
appendix, we have shown the mass spectrum of the singlet sector after minimizing the 
singlet potential. It turns out that the masses of axino and saxion partners in X, Y singlets 
are heavier than the gravitino mass. The details for the axion property are explained in 
Appendix B. Even with an explicit PQ-breaking term and after the i?-symmetry breakdown 
for nonzero singlet VEVs, the PQ symmetry breaking is small enough for maintaining the 
axion solution to the strong CP problem. The details on this aspect are shown in Appendix 
C. 



2.5 Neutrino masses 

In this section, we introduce right-handed neutrinos which can have Majorana or Dirac 
masses. It further constrains the allowed -R-charges of the MSSM sector, consequently 
determining the scalar soft masses via the U(1)r mediation. We first consider the Majorana 
neutrino case with intermediate or TeV-scale Majorana mass for the RH neutrino. Then, 
we go on to discuss the Dirac neutrino case with vanishing Majorana mass for the RH 
neutrino. 



2.5.1 Majorana neutrino case 

We consider the neutrino mass term in the superpotential with right-handed neutrino N 
as follows, 

W v = \ U LH U N + -^X n NN. (9) 



The standard high-scale see-saw mechanism is applied for n = 1 while the TeV-scale see- 
saw mechanism must be used for n = 2. Then, for the .R-charge of the Y singlet, ry = —3, 
we obtain the .R-charges of the other singlets as rx = § and = 1 — Then, the 
.R-charge of the squark doublet, that determines all the other .R-charges through Eq. (jHJ), 
is determined to be 

29 In 

« = -27-m- (10) 

Thus, we obtain the doublet squark .R-charge to be q = — 1| for n — 1 and q — — | for 
n = 2. Note that the PQ-charges can be assigned in an appropriate way and thus the axion 
solution to the strong CP problem persists even after introducing right-handed neutrinos. 
For n = 2, the i?-charges and PQ-charges are shown in Table 1. 



4 It is not possible to write the self-interaction for the singlet in the superpotential to break PQ symmetry 
explicitly. 

5 If the \i term comes from a renormalizable singlet interaction, a necessary small singlet VEV would 
lead to a dangerous axion due to low PQ symmetry breaking scale. 
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Table 1: R-charges and PQ-charges for the Majorana neutrino case with n = 2 

The PQ symmetry is nothing but a global /^-symmetry with new .R-charges fj given 
by the shifted ones from the local -R-charges, f i = r i J r qi with ^ being the PQ charges. As 
shown in the Appendix C, the PQ symmetry is broken explicitly by higher order U{1)r- 
invariant terms in the superpotential. The same is true of the Dirac neutrino case. 



2.5.2 Dirac neutrino case 

We note that in the absence of the Majorana neutrino mass term, it is possible in our 
framework to realize the tiny Dirac neutrino Yukawa coupling for neutrino masses in the 
superpotential. 

One possibility (Type I) is to take the following Dirac neutrino Yukawa coupling, 



W„ = ^XYLH U N. 



(11) 



(12) 



Since the -R-charge of the right-handed (RH) sneutrino becomes 

45 194 

the -R-charges of all fields except X, Y are determined in terms of the doublet squark R- 
charge. We note that the doublet squark .R-charge is not determined unlike the Majorana 
neutrino case. For this type of Dirac neutrino mass, taking into account the condition on 
q for getting positive soft squared masses of squarks, sleptons and RH sneutrino, we give 
an example with rational -R-charges in Table 2. 
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Table 2: R-charges and PQ-charges for the Dirac neutrino case 
The other possibility (Type II) is to take the following superpotential; 



Then, the -R-charge of the RH sneutrino is 



45 32 



(13) 



(14) 



We note that the Y 2 LH U N coupling would lead to the .R-charge of the RH sneutrino as 
r^r = y(? + ^p, giving rise to a tachyonic RH sneutrino for the allowed range of q. 

In either case without tachyonic RH sneutrino, a necessary tiny Dirac Yukawa coupling 
can be generated when the singlets get intermediate-scale VEVs as shown in Appendix A. 
For instance, in the former case, we obtain the neutrino mass as 

m„ = y u v sin ft ~ 0.01 eV , (15) 

where y u = jfr(XY). Thus, plugging the singlet VEVs (I A. 71) in the above, we require 

^7 ~ 10 4 for ~ 10~ 16 . In Type II case, we would need ^ ~ 10 3 for the same gravitino 
mass. 



3 Initial soft mass parameters 

Let us first summarize the [/(l)R-mediated SUSY breaking in the MSSM sector. Then, 
we also determine the soft mass parameters for the singlet sector, that is responsible for 
generating the \x term as well as the neutrino masses. 

For a scalar field with .R-charge r», the U(1) R D-term determines the scalar soft mass [9] 
as m\ = —riTriy 2 . Thus, from the -R-charges of the MSSM fields in Eq. fl5]), the scalar soft 
masses are given by 

m 2 ~ = -q m l /2j m 2 = (3$+ —^ml /2 , 
2 /3_ 26\ 2 2 ~ 18\ 2 2 ^31_ 46\ 2 

m e = (^ + — ) m 3/2^ m u=-(—(l+ — ) m 3/2^ m ^=(—(l+ — )m 3 



m 



' ^'"3/2) ni u - ~\y YJ /2 ' J ~ vT YJ /2 ' 

/24 60\ 2 2 /24 4\ 2 

h d = ~{- ( l + -) m V^ rn~ hu = {—q + -)m 3/2 . (16) 



Note that all squarks and leptons squared masses are positive when the doublet squark R- 
charge lies in the range — ||<g<— j|. In this i?-charge range, the soft mass squareds of 
the scalar Higgs doublets are negative. The corresponding trilinear soft terms for Yukawa 
couplings are universal as 

A ijk = -2m 3 /2 (17) 

for all i, j, k. 

From the modified gauge kinetic term, Eq. (JTj), in the presence of the nonzero F-term 
of the T-modulus, the gaugino masses for the SM gauge group are given by the £7(1)r-SM 
mixed anomalies and are universal at the GUT scale: 

M a = k a g 2 a F T ~ -—^m 3/2 , a = 1, 2, 3, (18) 

where the used relations are \Qn 2 g R k a g 2 = — 9(7q UT ~ — | and F T ~ 2m 3 / 2 - The gaugino 
masses can be larger or smaller compared to the gravitino mass depending on the ^(1)_r 
gauge coupling (e.g. \M a \ > m 3/2 for g R < j^). In particular, for a small U(1) R gauge 
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coupling as required for a large gaugino mass, the Green-Schwarz terms with large k a give 
a large negative contribution to the SM gauge kinetic terms. In this case, to get the unified 
value of the gauge couplings, .(/gut — |, we need to cancel the large contribution of the 
Green-Schwarz term by considering small tree-level SM gauge couplings. Henceforth we 
treat the universal gaugino mass My 2 to be a free parameter. We will see that My 2 3> m 3 / 2 
is required for a proper electroweak symmetry breaking. 

3.1 Majorana neutrino case 

For the Majorana neutrino case, the soft mass terms for singlets are 

£ so ft D -m 2 x \X\ 2 - m 2 Y \Y\ 2 - m 2 N \N\ 2 (19) 



-A h Y 2 H u H d - Tr A K YX i - \ U A U LH U N - —^ Z jA N X n NN + c.c. 



M P M P 2M r p 

The soft mass parameters for the singlet sector are determined in the U(1)r mediation as 
follows: 

m x = ~o m 3/2, «y = 3m 3/2 , = ( - 1 + -^) m 3/2 • ( 20 ) 



2 2 J/z: " ~iV y 6 

The A terms in the neutrino sector also follow the relation A ~ — 2m 3/ / 2 as shown below; 



/ h \ F c / F s F T 

Ah = - Fa ' ln (cTM:)^ +4 (67-ir 



-2m 3/2 , 



3/2, 



6s 3t 



-2m 3 / 2 , 



where ^ ~ | m 3/2 an d ^ — 2m 3 / 2 . Here, Yi's are defined from the expansion 
of the superconformal factor, Q = —3e~ K / 3 : Q ~ — 3e~ K °^ 3 + Y%Q\Qi where Qi are all 

the brane-localized chiral superfields, Kq is independent of the brane fields, and Yi = 

( \ 1/3 / \ -2/3 

( 2 (S + S^) ) 1 1 (T + T^) — Q'^Q' — (p^(p ) . In the presence of nonzero singlet VEVs, 
we obtain /i, B[i terms as follows, 

" = w P ^ (21 > 

Bp. = A h p + ^(YX- 3 ) =„[A h + ^-{Y- 1 X")). (22) 
After the singlet VEVs (1A.7j) are inserted in the above, we find that 

fi ~ 0.0272 -m 3/2 , (23) 
~ 7.49 m 3/2 . (24) 
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Moreover, the RH neutrino masses are also determined as follows, 



M N 
B N M N 



A N M N + ^^(Y*X**X^) 
Mn ( An + T^ {Y * X * 2X ' 1) )- 



(25) 



(26) 



In the n = 2 case, from Eq. ( TT6|) with i?-charges given in Table 1, the MSSM scalar soft 
masses at the GUT scale are determined as follows, 



m- 



2 ^ 2 2 

m i = 3^3/2, m i 



2 2^2 
m d = m e = 3 m 3/2, 



m 



m 



Ha 



-4^3/2- 



In this case, from Eqs. ( 125]) and ( 126]) . the RH neutrino masses are 

M N ~ 0.849 — m 3/2 , 
-Bjv ^ -1.09m 3/2 . 

Then, the mass eigenvalues of the RH sneutrino are 

m % = Ml + m 2 N ± \Bn\M? 



N± 



m 



3/2 



2 



0.849^) + -±( 0.849^ 
k I 3 \ K 



. 5 2 

> ^ m 3/2- 



(27) 



(28) 
(29) 



(30) 



3.2 Dirac neutrino case 

In the Dirac neutrino case with XYLH U N, the soft mass terms for singlets are 



£soft ^ 



-m x \X\ 



m z Y \Y\ 



m 



N 



h 



A h Y H u H r i - 



K 



-A K YX - 



K_ 

Ml 



A U XYLH U N + c.c. 



(31) 



The scalar soft masses for the X, Y singlets and the trilinear couplings corresponding to the 
(i term and the singlet interaction are the same as in the Majorana neutrino case, so, after 
X, Y singlets get VEVs, the \x and B\i terms are given by Eqs. (IE]) and ( 1241 . respectively. 

As the -R-charges of all fields are determined in terms of the squark i?-charge, so are 
the scalar soft masses in the U(1)r mediation. Then, the RH sneutrino scalar soft mass is 

given by m? N = — ^y<? + ^f^j m \/2- For the doublet squark i?-charge, — |f < q < — not 
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only all the squarks and sleptons but also the RH sneutrino have positive scalar squared 
soft massed as 



< m 2 N < 0.301 ml /2 . (32) 

We note that, in this region of \q\, the relic density coming from neutralino as the LSP 
tends to be too large [9j. So, it is natural to take the RH sneutrino or gravitino as a dark 
matter candidate. The trilinear coupling for the Dirac neutrino Yukawa coupling is given 
by 



2 ^ + ^-u^- 2m «'< 



following again the relation ~ — 2m 3 / 2 . 



4 Low energy spectrum and dark matter 

In this section, we consider the constraints on the SUSY spectrum coming from the EWSB 
conditions. Even after the i?-symmetry breakdown, the i?-parity is a good symmetry at 
the perturbative level as the i?-parity violating terms appear at sufficiently higher orders 
as shown in Appendix D. Depending on the nature of neutrino masses, we take either 
gravitino or RH sneutrino to be a dark matter candidate. We discuss the dark matter relic 
density and the BBN constraints on a late decaying NLSP in either dark matter scenario. 

4.1 The EWSB condition and the SUSY spectrum 

The Higgs mass terms contributing to the Higgs potential are given by 

V h , mass = (|/i| 2 + m 2 Hu )\H u \ 2 + (H 2 + m 2 Hd )\H d \ 2 + (B^H u H d + c.c.) (33) 

In order to achieve electroweak symmetry breaking, the following conditions at the weak 
scale must be fulfilled: 

|5H 2 >(H 2 + ^J(H 2 + m^), (34) 

2\n\ 2 + m 2 Hu +m 2 Hd -2\B f i\>0. (35) 
Then, the minimization conditions for the Higgs potential impose the following conditions; 

sin(2/3) = — 2 2 ^ - - , (36) 

I ,2 _ m H d ~ m 2 Hu tan 2 g M 2 

^ ~ tan 2 p - 1 T" - [ 7) 



6 For the Dirac Yukawa coupling X 2 LH U N, the RH sneutrino mass is given by 2.2Ani^ /2 < m 2 N < 
4.97m2 /2 . 
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The above EWSB conditions require that the /x-term and Higgs scalar soft masses must be 
large at the EWSB scale to be compatible with the large B-term as compared to gravitino 
or scalar soft masses at the GUT scale in Eq. (|24"I) . Therefore, the necessary large loop 
corrections to the Higgs scalar soft masses can be obtained for the gaugino mass which is 
much larger than the gravitino mass. 

When M1/2 3> m 3/2, after RGE running from GUT scale to EWSB scale, the soft 
terms at the EWSB scale becomes of the order the gaugino mass at GUT scale while B 
and /i-terms do not change much. For tan/3 > 1, considering |/i| 2 ~ — ~ ^1/2 from 
Eq. (137|) . we find roughly 



B ~ 7.5m 3/2 ~ J^L „ Ml/1. (38) 
tan fj tan p 



at the EWSB scale. This is the common feature of this U(1)r gauged model. For the 
correct magnitude for /i-term we need a large ratio between h and k as Ii/k ~ 10 3 from 
Eq. ©. 

Since gaugino mass is much larger than scalar soft masses at GUT scale, at low energy, 
all the masses of the SUSY particles are of the order of the gaugino mass and the lightest 
ordinary supersymmetric particle becomes lighter stau f\. Therefore the possible Dark 
Matter (DM) candidate must be gravitino or sneutrino which is outside of the MSSM sector. 
We note that the axino and saxion partners of the X, Y singlets are heavier than gravitino 
or sneutrino so they cannot be LSP. In the following sections we consider the corresponding 
DM candidate for each model introduced in the previous section. 

4.2 Majorana neutrino case 

First, we consider the Majorana neutrino with n = 1 in Eq. ((9]), the mass of the scalar RH 
neutrino mass is M N ~ \n(X) ~ 10 10_12 GeV. Thus, gravitino, as the LSP, is the only 
dark matter candidate while stau is the NLSP. In this case, there are two sources for the 
relic density of gravitino DM: non-thermal production from the decay of stau NLSP and 
thermal production from the thermal scattering after reheating. We will not consider the 
thermal production which depends on the reheating temperature. 

When NLSP decays after freezeout, the non-thermal production of LSP dark matter is 
determined by 

n DM h* = ^Lfi NLSP /i 2 . (39) 

^NLSP 

The abundance of stau from the thermal freeze-out is 



For gravitino LSP and stau NLSP which decay via fj — > r + G, from Eqs. fl39|) and ( ]4"0]) . 
the non-thermal production of gravitino is 
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For a correct relic density, we need a stau mass of TeV scale which is consistent with the 
EWSB conditions in our scenario. 

However the decay products of a long-lived decaying particle can be problematic with 
respect to the standard BBN. When a long-lived particle is negatively charged, it can make 
a bound state with nuclei and even worse the situation (CBBN) [25]. To avoid the CBBN 
constraint, the lifetime of stau must be less than 5 x 10 3 sec or Y fl = n fl /s < 1CT 15 for 
longer lifetime [2&H30] . Such a small abundance of stau requires unusual situations [3TH31] . 
As stau decays dominantly to gravitino and tau lepton, the lifetime of stau is given by 



r [f x ->• G + r ) ~ 48ttMJ> 



2 

2 m 3/2 



^ 1 . 8xl0 3 8eC (^-) 2 ('— V 

VIOOGeV/ V rrih J 

The CBBN constraint is automatically satisfied for > 2 TeV in the region of producing 
the right relic density of gravitino from stau decay according to Eq. (HTJ) . For example, for 
tan/3 = 10, our model with correct EWSB predicts m 3 / 2 — 200 GeV and m fl = 2.6 TeV, 
which leads to the lifetime of stau around 100 sec. 

In the case of Majorana neutrino with n = 2 in Eq. (Q, the mass of RH sneutrino 
is determined by Eq. (|30|) . If we consider the small magnitude of k ~ O(10~ 3 ) with 
An ~ (9(1); then m^ ± ~ (0.85Atv/k)?77 3 /2 ~ 10 3 m 3 / 2 and becomes much heavier than the 
other SUSY particles as well as gravitino. This gives gravitino DM as in the n = 1 case. 
For A at ~ ft, we may obtain ~ 777,3/2 so two scenarios are possible: RH sneutrino is 
LSP and gravitino is NLSP and vice versa. In both cases, however, the thermal production 
of the RH sneutrino LSP would highly overdose the Universe [23| , and thus it is excluded. 

In Figure [H we show the gaugino mass M1/2 at the GUT scale and the stau mass at the 
EWSB scale vs the gravitino mass after taking into account the correct EWSB. The green 
lines show the region where the gravitino non-thermal production from stau decay is within 
the range of cold dark matter from WMAP 7-year data [35], 0.105 < QoMh 2 < 0.119. The 
region right to the solid green line, where flDAih 2 > 0.119, is excluded. In Figure 2, we 
show the region of correct relic density in the plane of tan/3 and m 3 / 2 . 

Gravitino can be produced thermally, depending on the reheating temperature after 
inflation. In the region left to the green lines, the thermal production of gravitino is 
required for the correct relic density of gravitino. 



4.3 Dirac neutrino case 

First, we consider Dirac neutrino Type I with the superpotential Eq. (jTTj) . Since the 
renormalization group evolution of the (Dirac) RH sneutrino mass is negligible due to 
the smallness of Yukawa coupling, the RH sneutrino mass at EWSB is smaller than the 
gravitino mass from Eq. (13"2l . Therefore the RH sneutrino is the LSP Dark Matter, and 
gravitino is the NLSP. 
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Majorana neutiino with n=2 and tanfl=2 



m stau (EW) 



Majorana neutrino with n=2 and tanp=10 



m s „„(EW) 



100 
m 3/2 [GeV] 



m 3/2 [GeV] 



Figure 1: Plot of gaugino mass at GUT scale (upper blue line) and stau NLSP at EWSB 
(lower red line) for a given gravitino mass with correct EWSB condition in the case of Ma- 
jorana neutrino with n = 2. Between the vertical green lines, the non-thermally produced 
gravitino from stau decay satisfies the correct dark matter abundance by WMAP. 



Majorana neutiino with n=2 




50 100 150 200 250 

m 3/2 [GeV] 



Figure 2: The contour plot correct relic density of gravitino from non-thermal production 
in the plane of tan /3 and gravitino mass in the case of Majorana neutrino with n — 2. 



Let us note that the thermal production of the RH sneutrino can be obtained from the 
decay of supersymmetric particles but it is suppressed by the small Dirac neutrino Yukawa 
coupling y v ~ 1CT 13 . If there is no enhancement factor due to a small mass difference 
between left-handed and RH sneutrinos or degenerate neutrino masses (requiring large 
y u ), one typically gets the relic density of the RH sneutrino from the thermal production 
as flff ± h 2 < O(10~ 3 ) [23]. Thus, the main contribution to the RH sneutrino DM can come 
from the non-thermal production due to the decay of the NLSP. Using Eqs. f|59|) and ( HUj) 
for RH sneutrino LSP, the relic density of the RH sneutrino produced non-thermally from 
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Dirac neutrino type I with tan0=2 and qtilde=-1.48 (rN=-0.276) 



Dirac neutrino type I with tanp=10 and qtilde=-1.48 (rN=-0.276) 




m 3/2 [GeV] 



Figure 3: The same as Figured] but in the case of Dirac neutrino type I. 



Dirac neutrino type I 




50 100 150 200 250 

m y2 [GeV] 



Figure 4: The same as Figure [2] but in the case of Dirac neutrino type I. 



stau decay is 

" N ± h ~ m fl UTlh -^VlOOGevATeVy' [A6) 

which includes the RH sneutrino produced from stau decay to gravitino and gravitino decay 
to RH sneutrino. The region which gives the correct relic density for DM is shown in the 
Figure |3] for fixed tan/3 = 2, 10 and in the Figure H] on the tan/3 and 7713/2 plane. Here we 
used tn = —0.276 so that the mass of RH sneutrino is m^ ± = 0.52m 3 / 2 . With smaller 
the green lines in the Figures move to the right direction correspondingly. 

In the case of Dirac neutrino, the decay rate of stau to RH sneutrino and W gauge 
boson can be comparable to the one of stau to gravitino and tau lepton, causing a BBN 
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Dirac neutrino type I with qtilde=-1.48 (rN=-0.276) 



10° i|— 
Hi' f 
10" r 



i<r r 
10 1 r 
10° ' — 



m 3/2 [GeV] 



a io- 2 



Dirac neutrino type I with qtilde=-1.48 (rN=-0.276) 



m 3/2 [GeV] 



Figure 5: The lifetime (left) and the branching ratio of stau decay to RH sneutrino and W 
boson (right) for tan = 2, 10 in the case of Dirac neutrino type I. 



problem. The latter decay rate is given in Eq. (1421) and the former one is given by 

2 



r(fi ->• w~ + n±) 



sm 2 9 f ( m\ \ |/icot/3 - A* u \ 2 m 2 u 
32tt \ m2 u L ^ m \ ± J m fl v 2 



3.3 x IO" 26 GeV sin 2 9 f ( 1™) Y " (~^?) 

\m 0L J VlTeV/ \ 1 TeV J VO.OleV/ 



(44) 



2 



where use is made of v ~ 174GeV, is the neutrino mass, and 9? is the left-right 
mixing angle of stau, i.e. f\ = f R cos9 f + f L sin^. Here we used Eq. (fT5|) and A u ~ 
— 2m3/2 — 0.59Mi/2. Using this we show the plot of the lifetime of stau and the branching 
ratio of stau decay to RH sneutrino and W boson in Figure [5] for tan — 2, 10 respectively. 

As can be seen from Figure 5, for tan/3 = 10, one finds that the decay rate of stau to 
gravitino and tau lepton (|42|) is much larger than the one for stau to RH sneutrino and W 
gauge boson. Then, the lifetime of stau is determined by the decay rate to gravitino and 
tau lepton so it is less than about 100 sec. Therefore, in this case, the hadronic particles 
produced from W boson decay do not have a BBN problem [231124"] . For tan /3 smaller than 
10, the stau decay rate to RH sneutrino and W boson becomes sizable and thus it would 
cause the BBN problem. 

On the other hand, the late decay of gravitino to RH sneutrino and neutrino could cause 
a BBN problem. But, the BBN constraints on the late decaying gravitino may be avoided 
if the mass difference between gravitino and sneutrino is less than about 100 GeV [36J. 

For the Dirac neutrino Type II with the superpotential Eq. ([TBI , the mass range of RH 
sneutrino is 

1.49m 3 / 2 < m„ < 2.23m 3/2 . (45) 

In this case, RH sneutrino is NLSP and gravitino is LSP as a DM candidate. The correct 
relic density of gravitino can be obtained from the stau decay in some region of heavy stau 
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mass while the BBN constraint can be avoided in the same way as Dirac neutrino Type I 
model discussed above. 

In both Dirac neutrino cases, in the region where the non-thermal production is not 
enough, the thermal production of gravitino may give rise to the correct DM of gravitino (or 
RH sneutrino DM from gravitino decay) with appropriate reheating temperature around 
10 9_10 GeV (or scaled by ms/ 2 /m^). 

5 Conclusion 

We have shown that the gauged U{1)r symmetry naturally realizes the solution to the \x 
problem and accommodates the axion solution to the strong CP problem. The interplay 
between the higher dimensional interaction for singlets and the singlet soft masses coming 
from the U(1)r mediation gives rise to the stabilization of the singlets at an intermediate 
scale, consequently generating the small \x term from a higher dimensional interaction. 

The gauged U(1)r symmetry restricts the generated Bfj, term to be larger than the 
other scalar soft masses of order the gravitino mass, resulting in M 1 / 2 ^> m^n- Thus, we 
found that superpartner masses at the EWSB scale in the MSSM sector are much larger 
than the mass of gravitino or RH sneutrino. Therefore, only gravitino or RH sneutrino can 
be a natural dark matter candidate. 

Depending on whether the Majorana mass term for the RH neutrino exists, we con- 
sidered a different candidate for dark matter: gravitino for the Majorana neutrino case 
and RH sneutrino for the Dirac neutrino case. In both dark matter scenarios, the NLSP 
in the MSSM sector is stau. For the stau decaying after the freezeout, we showed that 
the correct relic density of dark matter can be generated by the non-thermal production 
mechanism through the stau decay. At the same time, the BBN constraints on such a 
long-lived charged particle can be evaded for the TeV-scale stau mass. Gravitino and RH 
sneutrino LSP could be also produced thermally, by the reheating after inflation and by 
the decay of the heavier superparticles in thermal bath, respectively. However, we have 
not considered the thermal production because the former case depends on the reheating 
temperature and the latter case depends on the (fine-tuned) enhancement factor of the 
decay rate. 
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Appendix A: Minimization of the scalar potential for 
singlets 
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In this section, in order to get nonzero singlet VEVs, we focus on the minimization of 
the scalar potential for the singlets, X and Y. 
The F-term potential for singlets X, Y is 



V F 



2h k 
■YH u H d + —-X 1 



M P M F 



Since the Higgs doublet get small VEVs compared to the singlet VEVs, we ignore the 
singlet coupling to the Higgs doublets in the F-term potential. As far as singlet F-terms are 
negligible and |Jf|, \ Y\ <C Mp, the U(1)r D-term contributes to the singlet scalar potential 
only through the soft mass terms, which are determined after the moduli stabilization. 
Then, adding the singlet soft mass terms to the F-term potential, we obtain the scalar 
potential for singlets as 

V(X,Y) = V F + m 2 x \X\ 2 + m 2 Y \Y\ 2 

- w P w° + S |y|2| - Y|4 + ' 4| - Y|2 + <]Yf + -k KYX * + C ' C - (A ' 2) 

Writing X = \X\e i9x , Y = \Y\e i6v , we find that for A K ~ -2m 3/2 < and k > 0, the 
trilinear terms stabilize one of linear combinations of angles at 36 x + @y = ^nn with integer 
n. Thus, the other combination of angles becomes a massless axion. After plugging the 
minimization condition for the angles into Eq. flA.2j) . the scalar potential becomes 

V(X, Y) = + g|y| W + m%\X? + - 2 -^\YUt (A.3) 

For the soft mass parameters given in the previous section, redefining the singlet fields as 

v 

V(x, y) = 3/2 P ( x 6 + 9y 2 x 4 - -x 2 + 3y 2 - Ayx 3 ) . (A.4) 
k V 3 / 

The extremum conditions for x and y are 

= 6x 5 + 36y 2 x 3 - —x - 12yx 2 , (A.5) 

3 

= 18yx 4 + 6y - 4x 3 . (A.6) 

Consequently, we find a minimum at x ~ 0.921 and y ~ 0.165 while x = y = is a saddle 
point. Then, the singlet VEVs are 



|2 



x 2 = ^ X [ f and y 2 = K ^ Y \, , we rewrite the scalar potential as 



\X\* 0.921 J m W M " t \ Y \* 0.165 J m ^ Mp . (A.7) 



K 



Expanding the singlets, X and Y, around the background VEVs, as X = (X) + + 

1^2) and Y = (Y) + ^75(^2 + ^2), we obtain the nonzero mass eigenvalues for singlets: for 
real bosons, 

M 2 h± = ~ (a + b ± v/(a - bf + 4c 2 ) (A.8) 
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with 



and 



ig<AT + %(Y?(X? - \ mil2 + ^-MY) (X), (A.9) 
^(X}* + 3ml /2 , (A.10) 



K - §f « 4 + if + ^ - (A - 12) 



for Weyl fermions, 



with 



X - [a! + V± - V) 2 + 4c/ 2 ) (A.13) 



«' - 5 <y)2<A ' >2+ 5 w4 ' (A - 14) 
= 5 <x>4 - (A - 15) 

c' = ^OW- (A.16) 

Another combination of the imaginary part is massless and it appears as a Goldstone boson 
for breaking the PQ symmetry. For the obtained singlet VEVs (1A.7jl . we can determine 
the mass eigenvalues: for the radial modes, which are almost mass eigenstates due to 
a small mixing, M| + = 9.Q7rriy 2 and M 2 _ = 8.39m3^ 2 ; for the massive angular mode, 
M 2 + = 12.2m\i 2 ; for Weyl fermions, M| + = 9.26m^ 2 and M|_ = 4.54m^ 2 . Here we note 
that the radial modes, h± are almost mass eigenstates /i2,i due to a small mixing. 



Appendix B: Axion for multiple scalar field VEVs 

We identify the axion when multiple scalar fields participate in PQ symmetry breaking. 

As shown in Appendix A, a linear combination of angles of singlet scalar fields, X 
and Y, in our model, i.e. 38x + Oy, is stabilized by the A-term for YX 3 term in the 
superpotential. From the PQ-charges of X and Y, this combination of angles does not 
transform under the U(1)pq. So, the orthogonal combination of angles plays a role for the 
QCD axion. 

From the kinetic term for X, Y, — (X) 2 (d ^6 x ) 2 — (Y) 2 (d^0Y) 2 , we find the canonical 
axion field as follows, 

a = jj((X)a x -3(Y)a Y ) (B.l) 
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where ax = jjk, a Y = Jyj an< ^ M = a/9(F) 2 + (X) 2 . In the presence of multiple scalars 
with VEV Vi and PQ-charge q\ the axion field is generalized [37J to a = -g J2i a % a%v i with 

Then, the axion coupling to the gluon field is given by the following effective Lagrangian, 

£a 99 = ^tr(G^) (B.2) 

where f a — ^ with A = J2i 1% being U(1)pq — SU(3)c — SU(3)c anomaly and U being the 
SU (3)c quadratic index of a fermion with PQ-charge q % . In our case, we obtain the anomaly 
as A = —3. Thus, the axion decay constant is given by |/ | = M/\A\ = |^/9(F) 2 + (X) 2 . 

Appendix C: PQ symmetry breaking terms and axion 
solution to strong CP problem 

We must also check other higher dimensional operators which are not PQ symmetric 
and thus can potentially spoil the property of the PQ symmetry for solving the strong CP 
problem. 

If there is a Planck-scale induced non-PQ symmetric term in the potential [38J: 

^ = ^0 2n+3 M+«*<n (ci) 

it gives additional contribution to the axion mass m 2 = |a||0| 2n+4 cos(5/(M|, r1 / 2 ) with 
5 being the phase of a. In order not to perturb the axion potential term from QCD 
instanton effect, this mass must be smaller than about 10~ 5 times the usual axion mass 
(m a ): m 2 < lCT 9 m 2 . 

First, we note that the PQ symmetry is an approximate global symmetry because 
it is broken explicitly by the Planck-scale suppressed £7(l)ij-invariant higher dimensional 
interactions, e.g. W = j^Y^X 12 . But, the leading term breaking the PQ symmetry 

while preserving the R symmetry is given by W = ■j^ n: (Wo) 4: Y 2 . Then, from the F-term 
potential for Y, we obtain the additional term for the axion potential as 

AV(a) = J ^(W rX* 3 Y + c. C . (C.2) 

Thus, by expanding the above potential around the axion minimum, we get the correction 
to the axion potential as 

/XV {a) ~ mla 2 - ^f a ml(tzn5)a (C.3) 
where m 2 = -^^{W ) A \X\ Z \Y\ cos 5 and e~ iS = ft. Then, from the bound (9) = & < 

yjaMp |0!| ' J a 

10~ 9 , we obtain 



m 2 tan 5 



ml + m 2 



< 1(T 9 (C.4) 
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I 2 1 

where m 2 n = Q S, D . Therefore, for Ira*! <C ra;?, this bound becomes < 10~ 9 . For 

nonzero singlet VEVs given in Eq. (IA.7I) . the axion mass correction is \m 2 \ ~ 0.3 x 

l a l("I7^r) m 3/2 ~ 10~ 26 eV 2 for m 3 / 2 ~ lOOGeV. Compared to the axion mass bound, 
m a = (0.6 x 10 7 GeV// a )eV > 0.6 x 10~ 5 eV for f a < 10 12 GeV, the axion mass correction 
is negligible. 

Appendix D: R-parity violating terms 

In this appendix, we consider the R-parity violating terms induced after the i?-symmetry 
breakdown. We focus on the Majorana neutrino case with n = 2. 

The effective R-parity violating terms are generated by the following higher dimensional 
interactions: 

1 YX 5 LQD, -LyX 5 LLE, -L(W ) 2 UDD, -±-W q Y 2 X 2 LH u . 



M e— Mf, ' My u; ' Ml 

However, after the superpotential and the singlets develop a nonzero VEV, the induced 
R-parity violating couplings are negligible. Therefore, it is possible to have a stable LSP. 
Moreover, we also note that the following terms are allowed: 

' (W ) 2 X 2 QQQL, -±-(W fYXUUDE. 



Thus, the B/L violating dim-5 operators are negligible so the proton stability is also jus- 
tified. 
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